let
  c = [-2., 1., 5.];
  G1 = [12. 13. 12.; 6. -3. -12.; -5. -5. 6.]';
  G2 = [3. 3. -1. 1.; -6. -6. -9. 19.; 10. -2. -2. -3.]';
  G = (G1,G2)
  h = ( [-12., -3., -2.],  [27., 0., 3., -42.]);
  sol = CVXOPT.socp(c, [], [], G, h)
end
